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Abstract 

We demonstrate nuclear double resonance for nanometer-scale volumes of spins where random 
fluctuations rather than Boltzmann polarization dominate. When the Hartmann-Hahn condition 
is met in a cross-polarization experiment, flip-flops occur between two species of spins and their 
fluctuations become coupled. We use magnetic resonance force microscopy to measure this effect 
between and ^'^C spins in ^'^C-enriched stearic acid. The development of a cross-polarization 
technique for statistical ensembles adds an important tool for generating chemical contrast in 
nanometer-scale magnetic resonance. 

PACS numbers: 76.70.-r, 05.40.-a, 76.60.-k. 76.60.Pc 
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The physics of microscopic spin ensembles can be distinctly different from that of macro- 
scopic ensembles. For example, in volumes of nuclear spins smaller than about (100 nm)^, 
random spin flips generate a fluctuating polarization that exceeds the typical thermal (or 
Boltzmann) polarization [1, 2, 3]. These spin fluctuations are a major source of dephas- 
ing in solid-state quantum systems [4, 5], and their control is an important prerequisite for 
nanometer-scale magnetic resonance imaging (MRI) and spectroscopy [6, 7, 8, 9, 10, 11]. Re- 
cent experiments using magnetic resonance force microscopy (MRFM) [12, 13] have extracted 
useful information from random polarization and harnessed it for nanometer-scale three- 
dimensional imaging [7] . One way to further improve nanometer-scale MRI is to combine its 
imaging capability with the chemical selectivity intrinsic to magnetic resonance. Techniques 
such as nuclear magnetic resonance (NMR) spectroscopy routinely exploit this feature for 
the structural analysis of molecules in millimeter-sized samples. At the nanometer-scale, the 
same information could be used to locally probe the chemical composition of materials, and 
to identify specific molecules in complex biological structures. 

Here we apply nuclear double resonance to achieve such a form of contrast using cross- 
polarization (CP) between statistically polarized and ^^C spins in ^^C-enriched stearic 
acid. CP is widely used in NMR for the signal enhancement of low-abundance and low-7 
nuclei and forms the basis for many advanced multidimensional spectroscopy techniques 
[14]. Indeed, CP has been demonstrated as an efficient chemical contrast mechanism for 
micrometer-scale one- dimensional MRFM imaging based on Boltzmann polarization [15, 16, 
17]. This spectroscopic method is not directly applicable to statistically polarized volumes 
of spins since at any given time the polarization has a random sign and magnitude, making 
measured signals intrinsically irreproducible. One way around this problem, as demonstrated 
here, is to observe the change in the correlation time of the fluctuations [3, 18]. 

CP relies on matching the rotating-frame Zeeman splittings of two different spin species 
(denoted as I and S) in order to promote cross-species spin flip-flops through the heteronu- 
clear dipolar coupling [14] . For this purpose, two strong rf fields are applied with frequencies 
near the Larmor resonance of the respective spins. According to the original work by Hart- 
mann and Hahn, efficient transfer is achieved when the rf field strengths are such that the 
respective Rabi frequencies have a similar magnitude [14, 19]. The Rabi frequency is deter- 
mined by the effective field in the rotating frame; for a spin /, the effective field is given 
by Beffi = Biix' + —l^vjz, where Bxj is the magnitude and Ai// = z// — is the reso- 
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nance offset of the rf field with a frequency z//. 77- is the nuclear gyromagnetic ratio of /, 
x' is a unit vector in the rotating frame, and z is a unit vector along the static field Bq 
(likewise for the S spin). The Hartmann-Hahn (HH) condition can then be expressed as 
7/ |Bcfr,i| = 75 |Bcfr,s|- The efficiency of CP depends on the angles of the effective fields with 
respect to the x' axis in the respective rotating frames, so that the most efficient transfer 
occurs when Az// = !\vs = and '^hBih = 'jcBic [15]. For samples with a high spin density 
and strong nuclear moments, the transfer process is typically very efficient, occurring on a 
characteristic time-scale set by the dipolar coupling frequency between spins. In stearic acid 
it is estimated at about 16 /is [20, 21]. 

For Boltzmann polarizations, a spin temperature description is commonly used in which 
the two spin polarizations are viewed as thermal ensembles [22]. This model successfully 
explains the dynamics of CP at the macroscopic level as well as in micrometer-scale MRFM 
[16]. When the HH condition is met, the two spin ensembles come into thermal contact and 
their temperatures equilibrate. In this way, one thermal ensemble can be used to enhance 
or deplete the polarization of the other ensemble. 

For statistically polarized nuclear spins, the same exchange of polarization occurs, except 
that the mean polarization of both spin ensembles is zero. When there is no HH contact, 
the polarization of an ensemble in the rotating frame fluctuates about zero with a variance 
that is given by the number of spins in that ensemble. The fluctuations in each ensemble 
occur on a time-scale determined by that ensemble's rotating frame relaxation rate. When 
the HH condition is met, double resonance allows rapid flip-flop processes that exchange 
polarization between the / and S spin ensembles, while simultaneously conserving the total 
polarization. As a result, fluctuations can occur on a much faster time-scale. 

To obtain a more quantitative picture, we consider the case for a spin-1/2 system where 
the polarizations of the ensembles are represented by rij and ns, defined as the difference 
between spin-up and spin-down populations for the I and S spins respectively. The dynamics 
of the HH transfer can then be described by a set of detailed balance equations that use 
spin population difference instead of spin temperature (see Ref. [19], Eq. (8)), 

^ = -kjnr - ^ {Nsnj - NjUs) + r,(t) + rjs{t), (1) 

= -ksns + ^ {Nsni - Nms) + rs{t) - ris{t), (2) 

where N is the total number of spins in the ensemble, Nj and Ns are the number of / and 



S spins respectively (A^ = Nj + Ns), kj and ks are the spin- lattice relaxation rates, and kjs 
is the average rate of exchange between I and S spins. 

The original work by Hartmann and Hahn (Ref. [19], Eq. (8)) does not account for 
statistical spin fluctuations. To describe the random excitations that lead to statistical 
polarization, we therefore introduce three stochastic functions rs{t) and rjs{t), equiv- 
alent to the number of /-spin flips, S'-spin flips, and cross-species spin flip-flops per unit 
time, respectively. We assume white spectral densities for these functions and further re- 
quire that on sufficiently long time-scales, the variance of the fluctuating polarizations is 
equal to the number of spins in the ensemble [1, 11]. The resulting (double-sided) spectral 
densities are Srj = 2kjNj and Sr^ = '^^ksNs- The spectral density Srjg = '^kjsNjNs/N is 
similarly obtained [23]. 

The coupled differential equations (1) and (2) can be solved in frequency space, yielding 
expressions for the spectral densities of rij and ns- One can then calculate experimental 
parameters, such as the variance or the characteristic time-scale of the spin fluctuations 
in the presence of CP. We concentrate our analysis on the spectral density of the J-spin 
fluctuations: 

^ m [{kj + kjs^) UJ' + {ks + kjs^) [kjks + kjs {kj^ + ks^)]] 

^4 + [^2 + ^2 + ^2^ + 2kjs {kj^ + ks^)] UJ^ + [kiks + kis (fc/f + ks^)f' 

Evaluation of (3) allows us to determine the behavior of the J-spin fluctuations for arbitrary 
values of /c/, ks and kjs- We consider two particularly relevant cases. As expected, in the 
regime of negligible CP where kis <^ /c/, ^s, the /-spin fluctuations occur on a time-scale 
Tj = kj^ with a variance equal to Nj: lim S", ^ ^ 



.^nj -. , 9 2' 

kjs/kj-,0 1 + UJ^Tf 

In the regime of strong CP where kjs 3> kj, ks, two time-scales emerge. Rapid polar- 
ization transfer between spin ensembles leads to fast spin fluctuations with a characteristic 
time Tis = kjg and a variance NjNs/N. The reduction of the variance compared to Nj 
occurs because the total polarization is conserved on time-scales short compared to r/ and 
Ts = kg^, thereby limiting the phase space of possible states. On top of the rapid exchange 
of polarization between spin ensembles, the polarization also fluctuates on a much slower 
time-scale Tavg = {kijf- + ksj^) ^ with a variance given by N] /N . Note that the sum of 
the rapid and slow variances is once again Ni. It is not surprising, then, that in this strong 
CP limit the /-spin spectral density approaches the sum of two spectral densities, each with 
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one of the two characteristic times and variances: 

(4) 



hm [Sni] = hm 



2tjsNjNs/N ^ 2Ta,gNf/N 



1 + LJ-^rfs 1 + uj^lg J 



Thus, the main signature of CP is the presence of spin fluctuations faster than tj. Simple 
coin-flipping simulations support these findings. 

We demonstrate nuclear CP between statistically polarized and ^^C spins in an ex- 
periment using a custom-built magnetic resonance force microscope [18]. For these species, 
where 'Jh/ic = 3.9772, the combination of the minimum Bu needed for adiabatic inver- 
sions and the maximum allowed current in our microfabricated rf field source prevent us 
from reaching the condition of most efficient CP, 'JhBih = IcBic- Nevertheless, we are 
able to observe significant CP at non-zero resonance offsets where the HH condition is met. 

The sample is a 10-/im-sized particle of stearic acid, C18H36O2, where >99% of the carbon 
is ^'^C. The particle is placed upon an ultrasensitive, single-crystal Si cantilever (120 //m long, 
3 yum wide, and 0.1 /im thick), with part of the particle sticking out beyond the end of the 
cantilever. The cantilever and particle are then briefly heated to a few °C below the melting 
point promoting adhesion to the cantilever. The cantilever with the sample attached is then 
mounted in a vacuum chamber (pressure < 1 x 10"^ torr) at the bottom of a cryostat, 
which is isolated from environmental vibrations. At the operating temperature of 4.2 K, the 
sample-loaded cantilever has a resonant frequency /o = I/Tq = 2.9 kHz, an intrinsic quality 
factor Qq = 44000, and a spring constant A; = 86 yuN/m. We actively damp the cantilever 
in order to give it a fast response time of ~ 25 ms. An FeCo nanomagnetic tip is used to 
produce the large (~ 10^ T/m) spatial field gradient required for generating magnetic forces 
of a few attonewtons between the spins in the sample and the tip. A microwire underneath 
the tip generates an rf field of a few millitesla that induces magnetic resonance in the sample 
[18]. The MRFM measurement is carried out on a stearic acid sample positioned 100 nm 
above the magnetic tip in an externally applied magnetic field |Be2;t| = 2.64 T. 

We measure the spin polarization by periodically inverting the nuclei of choice using 
adiabatic rapid passages [11, 18]. We operate in a fixed magnetic field Bq = ^ext + ^up 
(I Bo I = 2.72 T), where is the field produced by the magnetic tip. We periodically 
sweep the frequency uj of a transverse rf magnetic field Bi/ through the Larmor resonance 
condition, = so as to induce adiabatic inversions of the nuclear spin polarization. 

In the presence of the magnetic tip, periodic inversions of the spin polarization generate an 



alternating force that drives the mechanical resonance of the cantilever. The amplitude of 
cantilever oscillation, which we measure using optical interferometry, is then proportional 
to the /-spin polarization. As shown in Fig. 1, we can measure either the or the ^^C 
statistical polarization by adjusting the rf center frequency of the adiabatic passages. 

In order to observe polarization transfer, we perform adiabatic passages on one isotope — 
the "observed" or J-spin isotope. This measurement produces a signal that is proportional 
to the fluctuating /-spin polarization. Simultaneously, we address the other isotope — the 
"unobserved" or S'-spin isotope — with cw radiation at or near its Larmor frequency. This 
continuously applied resonant Bis is constant in the rotating frame and remains locked to 
a statistical polarization of S spins with a correlation time ts- The basic scheme is shown 
in Fig. 2. 

When the HH condition is fulfilled, polarization transfer occurs between the statistically 
polarized "observed" and "unobserved" spin ensembles. Measurements are shown in Fig. 3 
for both and ^^C in the role of the "observed" isotope. In both experiments we increment 
the resonance offset Aus of the S spins while detecting the / spins with frequency sweeps 
centered on Auj = 0. We record both the signal variance and the correlation time of the 
/-spin fluctuations. Dips in the correlation time appear at those frequencies where Az/g 
satisfies the HH condition. Since the signal is recorded in a narrow band (~ 17 Hz) around 
the cantilever resonance, a reduction in the correlation time due to CP also gives rise to a 
reduction in the observed force signal. 

Fig. 3 shows two distinct CP regimes. In the first regime, shown in Fig. 3(a), we observe 
the ^H spins, using 'JhBih = 280 kHz and 7c/?ic = 29 kHz, so that •yiBu > •ysBis- 
Again, as we increment Auq, 7c|Befjc| comes to a minimum at resonance according 
to 7c|Befr,c| = {icBicY + {'^t^AvcY. For two bands of Avc symmetric about zero, 
7c;|Befr,c| intersects the trajectory of 7_f/|Bcfr,H| thereby fulfilling the HH condition [16]. The 
resulting CP produces the double-dip structure shown in Fig. 3(a). The most efficient CP 
and therefore the most significant reduction in the ^H correlation time {i.e. the deepest 
part of the dips) occurs for a HH match at the vertex of the 7H|BefT,H| hyperbola shown 
in Fig. 3(a). At these intersections the slopes of 7c|Bcfr,c| and 7H|Befr,H| match, producing 
the longest possible HH contact. In addition, CP is most efficient when Avc and Av^ are 
smallest resulting in small angles of the effective field. 

In the second regime, shown in Fig. 3(b), we observe the ^^C spins using 'JhBih = 120 
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kHz and 7c-Bic = 62 kHz, so that 'yiBu < ^sBis- As we increment Auh, 7_f/|Bcff,H| comes 



there is one band of Auh symmetric about zero for which 7//|Bcfr,H| intersects the trajectory 
of 7c|Bcfr,c|- The resuhing CP produces the single-dip structure of Fig. 3(b). For the same 
reasons which apply in the first regime, the most efficient CP occurs for Auh = where the 
HH match is closest to the vertex of the 7c|Bcfr,c| hyperbola. The same single- and double- 
dip behavior shown in Fig. 3 appears in Fig. 2 of the original paper on CP by Hartmann 
and Hahn [19]. 

We have measured a series of spectra of the double-dip type for different 'JhBih and 
IcBic- Typical results are shown in Fig. 4(a). Using simple Lorentzian fits to determine the 
position of the dips, we extract the splitting for all spectra. From the condition 7/ |Bcfr,i| = 
7s'|Bcfr,s| and as pointed out by Hartmann and Hahn (Ref. [19], Eq. (75)), the dips in 
Fig. 4(a) should appear at 



so that the splitting is 2Avc,hh- The experimental data in Fig. 4(b) and (c) agree within the 
error with the theoretical curve representing (5) without any adjustable parameters. Note 
that the magnitudes of Bin and Bic are calibrated by independent nutation experiments 
on each species [18]. Given this agreement, we confirm that our double-resonance features 
result from the HH effect. 

While in our experiments it was possible to observe the ^^C signal directly, double- 
resonance detection is particularly useful when one isotope with spin S has a weak resonance 
signal, either because of a small 75 or a low abundance. The presence of the S spins 
can then be detected via the stronger signal of the / spins. The ability to perform CP 
in statistical ensembles provides the possibility of new contrast mechanisms for nanoscale 
MRI applications. For example, organic material with many proximate ^^C and ^H atoms 
could be distinguished from interstitial water molecules, which possess only ^H. Statistical 
double resonance could be combined with advanced spectroscopy techniques for chemical 
characterization of materials at the nanometer-scale. Ultimately, such techniques could 
be applied to discern individual protein components in complex nanometer-scale biological 
structures. 

We acknowledge support from the NSF-funded Center for Probing the Nanoscale (CPN) 



to a minimum at resonance according to 7// 1 Befr,H 



\/{'yHBiHy + {2tiAvhY- In this case. 
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at Stanford University. 
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FIG. 1: MRFM signal as a function of rf frequency at .Bq = 2.72 T for both the and ^^C 
isotopes in 99% ^'^C-enriched stearic acid. The two data sets were taken under shghtly different 
experimental conditions, making their relative amplitudes arbitrary. 
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FIG. 2: A schematic diagram showing the rf frequencies z// and us during an experiment. With 
vs constant, periodic sweeps of vi through resonance adiabaticahy invert the / spins. The effect 
of incrementing us is represented by the dotted arrow. 
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FIG. 3: CP resonances detected by MRFM in statistical polarizations plotted as a function of the 
(S-spin resonance offset. Below each plot, the corresponding trajectories of the effective field Rabi 
frequencies are shown. The effect of scanning is to vertically shift 75|Bcfr,s| as shown by the 
dotted arrow, (a) The "observed" spin / = and the "unobserved" spin S = ^^C with vpM = 300 
kHz. (b) The "observed" spin / = ^^C and the "unobserved" spin = -"^H with vpM = 125 kHz. 
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FIG. 4: (a) MRFM data (points) and fits (lines) as a function of the 5'-spin resonance offset for 
/ = and S = ^^C for different Bin, with Bic = 2.8 mT. (b) Ai^c,HH plotted as points for 
different Bih taken from fits in (a), (c) Avchh for different Bic from fits to a separate data set 
(not shown). Lines in (b) and (c) correspond to (5) with no free parameters. 
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